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Abstract 

Motivated by experiments in nanoscopic systems, we study a generalized Anderson, 
which consists of two spin degenerate doublets hybridized to a singlet by promotion 
of an electron to two conduction bands, as a function of the energy separation 5 
between both doublets. For 6=0 or very large, the model is equivalent to a one-level 
SU(N) Anderson model, with N=4 and 2 respectively. We study the evolution of the 
spectral density for both doublets (p\ a {uj) and p2 a (u})) and their width in the Kondo 
limit as 5 is varied, using the non-crossing approximation (NCA). As 5 increases, the 
peak at the Fermi energy in the spectral density (Kondo peak) splits and the density 
of the doublet of higher energy p2o-(w) shifts above the Ferrmi energy. The Kondo 
temperature Tk (determined by the half width at half maximum of the Kondo 
peak in density of the doublet of lower energy p% a (cj)) decreases dramatically. The 
variation of Tk with 5 is reproduced by a simple variational calculation. 
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1 Introduction 



The Kondo effect, found originally for systems with magnetic impurities in 
metals is now present in a variety of nanoscopic systems, including semicon- 
ducting quantum dots [1], magnetic adatoms on surfaces [1,2,3] and carbon 
nanotubes [4]. In the latter, in addition to the spin Kramers degeneracy, there 
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is in addition orbital degeneracy due to the "pseudospin" degree of freedom 
related with the particular band structure of graphene. This leads to an SU(4) 
Kondo effect which has been observed experimentally [5,6] and also discussed 
theoretically [7,8,9]. In particular, Lim et al. have studied the spectral density 
when the SU(4) symmetry is reduced to SU(2), mainly by a change in the 
tunneling matrix elements [7,9]. 

Our main motivation in the problem arises from interference phenomena in 
systems of quantum dots [10,11,12] or molecules [11,13,14]. For example de- 
pressions in the integrated conductance through a ring described by the Hub- 
bard or t — J model, pierced by an Aharonov-Bohm magnetic flux, related 
with spin-charge separation [10,11], are due to a partial destructive interfer- 
ence when the energy of two doublets cross [11]. Similar interference effects 
were predicted in molecular transistors [13,14]. The effective Hamiltonian near 
the crossing is discussed in the next section. 

To our knowlege, the conductance in these systems has so far been calculated 
using approximate expressions or a slave-boson formalism [12], which is valid 
only for very low temperatures and applied bias voltages. This work is a step 
towards a more quantitative theory to describe the transport through similar 
systems, treating the effective Hamiltonian within the non-crossing approxi- 
mation (NCA) [15,16]. Work is in progress to deal with the non-equilibrium 
situation, which is necessary within our formalism to calculate the current. 

In this paper, we report on our study of the spectral density of the model 
as a function of the splitting S between both doublets in the Kondo regime. 
We also calculate the dependence of the Kondo temperature Tk with 5 and 
analyze the validity of the Friedel sum rule [17]. 



2 Model 



We start from a model in which two doublets of an interacting system are hy- 
bridized with a singlet by promotion of an electron to two conducting leads. 
This is the low-energy effective Hamiltonian for several systems with par- 
tial destructive interference, such as Aharonov-Bohm rings [11] or aromatic 
molecules [13,14] connected to conducting leads. The Hamiltonian can be writ- 
ten as [11] 

ia uka 

+ £(K,R(o|c 1 , jfcCT + H.c.), (i) 

ivka 
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where the singlet |0) and the two doublets \icr) (i — 1,2; a =t or l) denote the 
localized states (representing for example the low-energy states of an isolated 
molecule), c vka create conduction states in the left (v = L) or right [y = R) 
lead, and V iu describe the four hopping elements between the two leads and 
both doublets, assumed independent of k. 

Changing the phase of the conduction states and the relative phase between 
both doublets, three among the four V iv can be made real and positive. The 
phase (p of the remaining hopping |T/|e 1 ^, depends on the particular system 
and its symmetry. For example in molecules with rotational symmetry = 
(K\ — K 2 )l, where I is the distance between the sites connected to the left and 
right leads, and Ki is the wave vector of the state \icr), which can be modified 
with an applied magnetic flux [11]. In absence of magnetic flux and when the 
relevant states have wave vector Ki = ±7r/2, as in rings with a number of 
atoms multiple of four, = n and there is complete destructive interference 
in transport [11,14]. 

In the following we will assume this case, with states 1 and 2 related by 
symmetry implying \Vi„\ = \V 2v \. We further assume symmetric leads, eik = 
e_Rfc = £fe, \Vil\ = \Vm\- Then, without loss of generality we can take Vu = 
V 1R = V 2L = V>0, V 2R = -V, and E 2 > E x . 

For these parameters, changing basis c\ ka = {c[ ka + c Rka )/V^, c\ ka = {c[ ka - 
c Rka )/y/2, the Hamiltonian takes the form of an SU(4) Anderson model with 
a "field" 5 = E 2 — E 1 and on-site hybridization V = \/2V 

ia ika 

+\/'^(|m)(0| Qfc(7 + H.c.). (2) 

ika 

Interchanging the doublet index (1 or 2) with the spin index a, one realizes that 
this model also describes transport trough a carbon nanotube with electron 
density depleated at two points (so as to created an SU(4) quantum dot in 
the middle) under a real applied magnetic field. 



3 Spectral density 

In Fig. 1 we present numerical results for the spectral density pi a of the SU(4) 
Anderson model (E\ = E 2 ) in the Kondo regime ep — ^> A, where the hy- 
bridization function A = n I]fc(V) 2 5(w — e^), assumed independent of energy. 
We set T = 2A = 1 as the unit of energy. We also assume a conduction band 
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Fig. 1. Spectral density as a function of energy for different temperatures. The 
inset shows a detail near the Fermi energy. Parameters are A = 0.5. D = 10, 
Ei = E 2 = -4. The lowest temperature is T = 10~ 3 = 0.0761^. 




Fig. 2. Spectral densities for levels 1 (full line) and 2 (dashed line) as a function of 
energy for different values of E 2 = E\ + 5 and T = 10 -3 . Other parameters as in 
Fig. 1. Dot-dashed line corresponds to p2 a for 5 = 0.015, 0.3 and 0.6. 



symmetric around the Fermi level ep = 0, of half width D = 10. 



The spectral density shows a broad charge transfer peak near E]_. For tem- 
peratures below a characteristic energy scale T K (defined below), p ia develops 
a narrow peak around the Fermi level. In contrast to the better known one- 
level SU(2) case, this peak is displaced towards positive energies and is much 
broader, as discussed in Section 5. 
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The evolution of the spectral densities as E 2 is displaced to larger energies, 
breaking the SU(4) symmetry is shown in Fig. 2. The peak near the Fermi 
energy of P2a{< jJ ) is displaced towards positive energies (near 5 — E 2 — Ex). In 
contrast, the corresponding peak in pi a (oj) narrows significantly and displaces 
towards the Fermi energy. This implies that the Kondo temperature TV defined 
as the half width at half maximum of this peak, also decreases strongly. The 
evolution of TV with 8 is discussed in the Section 5. 

In addition p lcr develops a broad peak near energy —5 which becomes visible 
when 5 becomes greater than T K . 



4 Friedel sum rule 



The Anderson model studied has a Fermi liquid ground state which satisfies 
well known relationships at zero temperature. One of them is the Friedel sum 
rule which relates the spectral density at the Fermi level for each "pseudospin" 
channel with the occupation of that channel [17]. For the simplest case of a 
constant density of conduction states, this rule reads 

pia{tF) = ^T-sin 2 (7rn i(T ), (3) 



where n ia = (\ia) (ia\) . 

This is an exact relationship for a Fermi liquid, which is not necessarily sat- 
isfied by approximations. In particular, it is known that at very low temper- 
atures, the NCA has a tendency to develop spurious spikes in Pi a {uS) at the 
Fermi energy, while thermodynamic properties, such as expectation values are 
accurately reproduced [15,16]. 

In Fig. 3, we compare both members of Eq. (3) for the lowest lying doublet, at 
a temperature T = 0.17V [18] low enough so that no further increase in pi a (ep) 
takes place as the temperature is lowered (according to physical expectations 
and Eq. (3)), but high enough to prevent the presence of spurious spikes. The 
disagreement lies below 20 %. The agreement improves as the parameters are 
moved deeper in the Kondo regime e F — E x 3> A. Thus, while the spectral 
density at zero temperature is not well represented by the NCA results at 
T = 0, one can take the values at T = 0.17V as a reasonable description of 
the correct T = ones. This statement is supported by the comparison of 
results obtained by NCA and numerical renormalization group for the one- 
level SU(2) case [19]. 
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Fig. 3. Squares: rescaled spectral density of the lowest lying level at the Fermi energy 
as a function of 5. Circles: corresponding (more accurate) result given by Eq. (3). 
Lines are guides to the eye. Parameters are T = O.ITk and the rest as in Fig. 1. 

5 Kondo temperature 



From the half width at half maximum of the peak nearest to the Fermi energy 
of the spectral density of the lowest level (p lo .(u;)), we have calculated the 
Kondo temperature of the system Tk for several values of 5. This requires to 
solve the self-consistent NCA equations up to low enough temperatures (about 
0.1 Tk as discussed above) so that the height of the peak does not increase 
significantly with further lowering of the temperature [18]. Fortunately, the 
result is not very sensitive to the ratio T /Tk- 

The results are shown in Fig. 4 and compared with Eq. (5) obtained from a 
variational calculation as explained below. We see that except for an overall 
multiplicative factor, the agreement is quite good, in spite of the fact that T K 
changes by nearly two orders of magnitude. 

To provide an independent estimate of T^, we have calculated the stabilization 
energy of the following variational wave function 

|V>) = a\s) +J2Miv)(0\c ik *)\s), (4) 

ika 



where \s) is the many-body singlet state with the filled Fermi sea of conduction 
electrons and the state |0) at the localized site, while a and {5$. are variational 
parameters. From the resulting optimized energy E, we can define the stabi- 
lization energy as T K = E\ — E. The Kondo energy scale defined in this way 
becomes 

T* K = {(D + 5)Dexp [nE 1 /(2A)] + 5 2 /a} 1 ' 2 - 5/2. (5) 
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Fig. 4. Squares: Kondo energy scale determined by the width of the peak in the 
spectral density near the Fermi energy as a function of the splitting 5. The temper- 
atures used were T = 10~ 3 , T = 10~ 4 , and T = 10 -5 depending on S Dashed line: 
corresponding variational result Eq. (5) multiplied by a factor 0.606. 



This expression interpolates between the SU(4) result (T£ = D exp [irE 1 / (4A)] 
for 5 = 0) and the SU(2) one for one doublet only {T* K = Dexp [ttEi/(2A)] 
for S — >• +oo). 



6 Summary and discussion 



We have studied an impurity Anderson model containing two doublets, which 
interpolates between the cases for one level with SU(4) and SU(2) symme- 
try, and is of interest for several nanoscopic systems, using the non-crossing 
approximation (NCA). 

We have shown that the NCA provides reasonable results for the equilibrium 
spectral density. The values of the spectral density for both doublets agree 
within 20% with the predictions of the Friedel sum rule, in spite of the fact 
that it is not expected to satisfy Fermi liquid relationships at zero temperature. 

In addition, the Kondo temperature scale Tk obtained from the width of 
the peak in spectral density near the Fermi energy agrees very well with the 
stabilization energy of a variational calculation, in spite of the change in sev- 
eral orders of magnitude of T% when the splitting between both doublets is 
changed. 

The approach seems promising for studying transport properties within a non- 
equilibrium formalism. Work in this direction is in progress. 
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